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\ Abstract. We present preliminary results from the Shellfiow program, 

■ an all-sky Tully-Fisher (TF) survey of 297 Sb— Sc with redshifts between 
^ ! 4000 and 7500 kms'^ The program was designed to ensure uniformity of 

J> ' the data between observing runs and between telescopes, thereby elimi- 

*^ , nating the possibility of a spurious bulk flow caused by data inhomogene- 

' ity. A simple bulk flow model suggests that the Shellfiow galaxies are 

. nearly at rest with respect to the cosmic microwave background (CMB). 

I Taken at face value, this result suggests that most of the ~ 600 kms"^ 

On . motion of the Local Group with respect to the CMB is due to mate- 

I rial within ~ 60/i~^Mpc — a striking confirmation of the gravitational 

. instability picture and the notion of homogeneity in the Universe. Dif- 

Qh| ferences between our Shellfiow "null" result and claims of larger bulk 

Q ■ motions based on analyses of the Marklll catalog are possibly due to 

|H \ small inhomogeneities between the sub-samples in Marklll. Our prelimi- 

c/3 . nary Shellfiow bulk velocity, V\)uik = 80 it 150 kms~^ in the CMB frame, 

' must still be refined with Monte-Carlo simulations and tidal field analysis. 
> ' 

. F~ , 1. Introduction 

It is of great cosmological importance to identify the volume of space, centered 
on the Local Group, which is at rest with respect to the CMB. Very large scale 
fluctuations are required to move ever larger volumes of space in the standard 
gravitational instability model. In standard Cold Dark Matter cosmologies, the 
volume of space bounded by the nearest superclusters (Great Attractor, Pisces- 
Perseus, Coma) is expected to define an inertial frame nearly at rest with respect 
to the CMB. The distribution of the matter within this volume should explain 
the ~ 600 km s~^ motion of the Local Group in the CMB frame. The detection of 
a large amplitude {Vhuik ~ 700 kms"^) on scales exceeding 8000 kms~^ by Lauer 
&; Postman (1994), followed by recent measurements with similar amplitude by 
Willick (these proceedings) and the SMAC team (Smith, these proceedings), 
have not only challenged the notion that the bulk flow on large scales is small, 
but have also pushed cosmological models to the breaking point. 
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Moreover, there are contradictory claims for the observed bulk flow within 
a sphere of 6000 kms~^, and whether it is generated by internal or external mass 
fluctuations. The most recent POTENT reconstruction of the Marklll velocities 
(Dekel et al. 1999) found a bulk velocity within 6000 kms"^ of 370 ± 125 km s"^ 
in the CMB frame towards Supergalactic {L,B) = (165°, —10°). Dekel et al. 
(1999) argue that this motion is generated by the external mass distribution on 
very large scales (see also Courteau et al. 1993). However, other investigations 
using nearly homogeneous samples of galaxies within and beyond ~ 60/i~^Mpc 
find motion consistent with the amplitude and direction of the CMB dipole (Gio- 
vanelli, these proceedings). New local peculiar velocity measurements based on 
the brightness of supernovae Type la (Riess, these proceedings) and surface 
brightness fiuctuations (Tonry and Dressier, these proceedings) are also consis- 
tent with little bulk motion at 60/i~^Mpc. This suggests that the reflex motion 
of the Local Group could be explained by material contained within 60 /i~^Mpc. 

The controversy over the observed bulk flow within 60/i~^Mpc stems, in 
large part, from our inability to combine the various galaxy distance samples 
used in flow studies into a single homogeneous catalog. None of the surveys 
within 60/i~^Mpc had sampled the entire sky uniformly]^. In an effort to over- 
come this situation, we devised a new study that would take advantage of North- 
ern and Southern NOAO facilities to map the motions of a shell of galaxies 
with signiflcant overlap between the hemispheres, and with existing peculiar 
velocity surveys. Our new ^^Shellflow" survey was designed to provide precise 
and uniform photometric and spectroscopic data over the whole sky, and thus 
to remove the uncertainties associated with matching heterogeneous data sets. 
With Shellflow, and further calibration of existing TF data sets, we expect a 
high-accuracy detection of the bulk flow amplitude and a better constrained 
characterization of the tidal fleld at 6000 kms~^. 

2. Observational Strategy 

The Shellflow sample is drawn from the Optical Redshift Survey sample of San- 
tiago et al. (1995; ORS). The ORS sample consists of all galaxies in the UGC, 
ESO, and ESGC Catalogs with tub < 14.5 and \b\ > 20°. We initially selected all 
non-interacting Sb and Sc galaxies in the ORS with redshiftsQ between 4000 and 
7500 kms~^, inclinations between 45° and 78°, and with Burstein-Heiles (1982) 
extinctions As <: 0™30. This produced a catalog of more than 300 objects. All 
galaxies were then visually inspected on the Digitized POSS plates; those with 
bright foreground stars and obvious disturbances were excluded, yielding a final 
sample of 297 Shellflow galaxies. No pruning was done of galaxies not matching 
idealized morphologies beyond the restriction on Hubble type and inclination. 



^Earlier attempts include David Schlegel and Josh Roth's thesis work. 

■^We actually define three subsamples complete in that range with diflFerent definitions of redshift: 
measured in the Local Group frame, the CMB frame, and after correction for peculiar velocities 
according to the IRAS model of Yahil et al. (1991) with /3 = 1; if we chose only one of them, 
the sample would decrease in size by 20%. 
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The data were collected between March 1996 and March 1998 at NOAO, 
using V and I-band photometry and Ha rotation curves. Data taking and reduc- 
tion techniques follow the basic guidelines of previous optical TF surveys {e.g. 
Courteau 1996, 1997; Schlegel 1995). The V and I-band imaging was obtained at 
the CTIO and KPNO 0.9m telescopes. The photometric calibration is based on 
the Kron-Cousins system; data taken on nights with standard star photometric 
scatter greater than 0'?020 were excluded. The Kron-Cousins system also al- 
lows direct matching with two largest I-band TF samples to date (Mathewson 
et ai. 1992, GiovancUi ct al. 1998). The Ha spectroscopy was obtained mostly 
in photometric conditions with the RC spectrographs at the CTIO and KPNO 
4m telescopes. Typical integrations were ~ 900s and ^ 1800s for imaging and 
spectroscopy respectively. 

Forty-one galaxies were imaged at both CTIO and KPNO, and we have 
repeat imaging from a given telescope for 106 galaxies. In addition, we observed 
27 galaxies spectroscopically from both CTIO and KPNO, and obtained dupli- 
cate spectra from a given telescope for 38 galaxies. The total magnitudes and 
rotational line widths reproduce to within 0"?06 and 3 kms~-^ (rms deviations) 
respectively, with no systematic effects seen between hemispheres or between 
runs. Moreover, all data reduction was done independently by Courteau and 
Willick using different software and methodology; the results between the two 
agree to within the errors quoted above. The high level of accuracy in our 
Shellfiow data meets our requirement for a measurement of a significant bulk 
flow result. 

3. Analysis and Results 

Following the approach of Lauer & Postman (1994), we use the sample itself 
to calibrate the distance indicator relation; this mitigates the need to tie the 
sample to external TF calibrators such as clusters (although it precludes mea- 
surement of a monopole term in the velocity field). We choose the "inverse" 
form of the TF relation, which makes Malmquist and selection bias eff^ects neg- 
ligible (cf., Strauss & Willick 1995). A maximum likelihood analysis for a pure 
Hubble expansion model (no bulk flow) allows us to constrain o"*"*, the intrinsic 
TF scatter, and 6v, the typical error in measuring raw rotation velocities. We 
hardwire Cp/iot, the total error associated with photometry (including inclination 
errors), to 0.15 mag, which is a small fraction of the error budget. Two other pa- 
rameters constrain the dependence of the TF relation on surface brightness and 
shape of the rotation curve. The best TF fit is achieved for Vrot measTired near 
1.7 disk scale lengths, similar to but somewhat smaller than the values found by 
Courteau (1997) and Willick (1999). A small surface brightness dependence is 
also found for the optimized TF relation, although this could be induced by the 
particular choice of luminosity profile fitting that we have adopted (see also the 
discussion by Willick [LPIOK], these proceedings; Courteau et al. 2000). 
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Shellflow TF Relation 




-23 -22 -21 -20 -19 

M + 0.35(^1-18.6) - 0.27(c-2.4) 



Figure 1. Tully-Fisher calibration of the Shellflow sample. This fig- 
ure includes all 297 Shellflow galaxies. The slope and scatter are com- 
parable to current modern I-band TF investigations. The absolute 
magnitudes are corrected for the small surface brightness and concen- 
tration index dependences of the TF relation. As mentioned in the 
text, the reality of these dependences is still under investigation by 
ourselves and other workers. The results of our flow fits are, how- 
ever, independent of whether or not these additional correlations are 
included. 
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For this analysis, we adopt a model of Hubble expansion plus a bulk 

flow to compute absolute from apparent magnitudes: 



d(Mpc) = ^(cz-VB-n) , 

where Hq = 100 kms~^ Mpc~^ and cz is the measured redshift in kms~^. We 
measure cz in either the CMB frame or the LG frame. (A different choice of 
Hubble constant would shift the zero point of the derived TF relation but would 
have no effect on our bulk flow fits.) 

The TF fit slope and scatter are comparable to recent I-band investigations 
{e.g. Giovanelli et al. 1998). The Shellflow TF relation is shown in Fig. 1. We find 
no dependence of the internal extinction correction on luminosity, and the TF 
parameters are virtually unchanged whether one adopts Burstein-Heiles (BH), 
or Schlegel, Finkbeiner, & Davis 1998 (SFD) Galactic extinctions. The opti- 
mization with respect to all TF parameters and bulk flow velocity components 
(Vx,Vjy,Vj;), yields the following fits: 

Table 1. Best-fit Velocity Components (kms~^) 

SFD extinctions 



Yz Yx Vy Frame 



47.5 -37.7 28.6 CMB 
-257.1 -98.9 566.1 LG 

BH extinctions 



Yz Yx Yy Frame 



58.8 -38.0 46.6 
-242.6 -99.2 584.5 



CMB 
LG 



The typical velocity errors, computed by holding the other parameters fixed, 
are AV^ = ±110 kms-^ AVj^ = ±90 kms-^ AV^ = ±70 kms-^ for an 
upper limit on the total amplitude equal to AV < 150 kms^^(lcr). Numerical 
simulations to assess the full covariance matrix of the errors will be presented 
elsewhere. 

Modelling the Shellflow galaxy peculiar velocities with a Hubble expansion 
and a bulk flow, one finds a low bulk motion of 80 ± 150 kms^^ with respect 
to the CMB. Correspondingly, the Shellflow volume moves at ~ 600 kms~^ in 
the frame of the LG^. This corresponds to the reflex motion of the Local Group 
motion with respect to the CMB. The directional components of the flow are 
poorly constrained since its amplitude is so low. 



■^See Courteau & van den Bergh (1999) for caveats regarding transformations to the Local Group 
reference frame. 
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Velocity Residuals, CMB Frame 




(point size « v ) filled circle: CMB apex 



Figure 2. Plot of peculiar velocities in the CMB frame. The point 
size is proportional to the velocity amplitude, as indicated by the key 
at the lower left. The circles and asterisks represent inflowing and 
outflowing objects, respectively. The two point types are well-mixed at 
all positions on the sky, indicating the absence of any significant bulk 
flow in the CMB frame. 
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Our result is in agreement with other newly-presented results at this work- 
shop, derived from surface brightness fluctuation techniques (Tonry h Dressier; 
Vfcu/fc = 289 lb 137 kms"-*^ at 3000 kms~^), nearby SNIa {Vbuik consistent with 
zero at 6000 kms~^(Riess, these proceedings), and previous TF analyses by 
Giovanelli and collaborators. Taken together these results suggest that most of 
the mass responsible for the motion of the LG lies within 6000 kms~^. Thanks 
to many new peculiar velocity surveys, the picture of cosmological bulk flows 
seems to imply "convergence" of the flow field within 6000 kms^^ to the CMB 
dipole value. However, this simple picture is challenged by a mix of low and 
high amplitude bulk flow measurements on larger scales. If the flow is nearly 
quiet at ~ 6000 — 7000 kms~^, one is hard-pressed to explain the reported high 
amplitude flows on larger scales. These issues are discussed elsewhere in these 
proceedings. 

4. The MarklV Catalog of Galaxy Peculiar Velocities 

Two of us (JW &; SC -I- collaborators) have recently combined the major distance- 
redshift surveys from both hemispheres (published before 1994) into a catalog 
of 3100 galaxies (Willick et al. 1997, the "Marklll Catalog"). However, fuU ho- 
mogenization at the 2 — 3% level, the minimum required for an accurate bulk 
flow detection at 6000 kms~^, was not achieved in Marklll due to irreducible 
uncertainties associated with matching disparate TF datasets. Furthermore, a 
revised calibration of the Marklll TF zero-points based on maximal agreement 
with the peculiar velocities predicted by the IRAS 1.2Jy redshift survey sug- 
gests a possible source of systematic error for the data sets which cover the PP 
cone (Davis, Nusser, & Willick 1996, WilHck & Strauss 1998 [VELMOD]). This 
uncertainty has not seriously affected mass density reconstructions (Dekel et al. 
1999) but it could lead to spurious estimates of the bulk flows on larger scales. 

Marklll was originally tied to the North/South cluster calibration of Han 
& Mould (1992). We now plan a revised MarklV calibration based on Shellflow. 
Because of the overlap with existing surveys at comparable depth, Shellflow will 
be of fundamental importance in tying these datasets together in a uniform way. 
With this more accurate global calibration, we expect that the MarklV catalog 
will be better sTiited than its predecessors for studying the velocity and density 
fields in the local universe. 

5. Summary 

• The Shellflow velocity components in the CMB frame are small. There is 
no significant bulk motion of the Shellflow sample relative to the CMB. 

• The Shellflow velocity components are large when the flt is done in the 
LG frame. This represents the reflex of the LG motion with respect to the 
CMB. 

• Results are insensitive to BH versus SFD Galactic extinctions. 
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• The Shcllflow bulk motion reported here has not been fully tested against 
Monte-Carlo simulations or using a complete tidal field analysis. A com- 
prehensive study is under way. 



Acknowledgments. We wish to thank various students and postdocs who 
have contributed to some of the Shellflow reductions. In Victoria, Shelly Pinder 
and Yong-John Sohn, and Josh Simon, Felicia Tarn, and Marcos Lopez-Caniego 
at Stanford. MAS acknowledges support from Research Corporation and NSF 
grant AST96- 16901. JAW acknowledges support from Research Corporation 
and NSF grant AST96-17188. 

References 

Burstein, D., & Heiles, C. 1982, AJ, 87, 1165 (BH) 
Courteau, S. 1996, ApJS, 103, 363 

Courteau, S. 1997, AJ, 114, 2402 

Courteau, S., Faber, S.M., Dressier, A., & Willick, J.A. 1993, ApJ, 412, L51 
Courteau, S., k van den Bergh, S. 1999, AJ, 118, 337 
Courteau, S. et al. 2000, in preparation 
Davis, M, Nusser, A., & Willick, J. A., 1996, ApJ, 473, 22 
Dekel, A., Eldar, A., Kolatt, T., Yahil, A., Willick, J. A., Faber, S. M., Courteau, 
S., & Burstein, D. 1999, ApJ, 522, 1 

Giovanelli, R., Haynes, M.P., Freudling, W., da Costa, L. N., Salzer, J.J., & 

Wegner, G. 1998, ApJ, 505, L91 
Han, M.-S., & Mould, J. R. 1992, ApJ, 396, 453 
Lauer, T. R., & Postman, M. 1994, ApJ, 425, 418 
Mathewson, D. S., Ford, V. L, & Buchhorn, M. 1992, ApJS, 81, 413 

Santiago, B. X., Strauss, M. A., Lahav, O., Davis, M., 8z Huchra, J. P. 1995, 
ApJ, 446, 457 

Schlegel, D. 1995, PhD. Thesis, University of California 

Schlegel, D., Finkbeiner, D.P., & Davis, M. 1998, ApJ, 500, 525 (SFD) 

Strauss, M.A., & Willick, J.A. 1995, Physics Reports, 261, 271 

Willick, J.A. 1999, ApJ, 516, 47 

Willick, J.A., & Strauss, M.A. 1998, ApJ, 507, 64 

Wilhck, J. A. et al. 1997, ApJS, 109, 333 

Yahil, A., Strauss, M. A., Davis, M., & Huchra, J. P. 1991, ApJ, 372, 380 



8 



